function y = shiftMDL(model,dr)

        for i=1:size(model,3),
            for j=1:size(model,2),
                model(1,j,i) = model(1,j,i) + dr(1);
                model(2,j,i) = model(2,j,i) + dr(2);
            end
        end

y = model;
end